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Abstract 

We present a new method, the analytical rebridging scheme, for Monte Carlo 
simulation of proline-containing, cyclic peptides. The cis/trans isomeriza- 
tion is accommodated by allowing for two states of the amide bond. We 
apply our method to five peptides that have previously been characterized by 
NMR methods. Our simulations achieve effective equilibration and agree well 
with experimental data in all cases. We discuss the importance of effective 
equilibration as well as the role of bond flexibility and solvent effects in the 
prediction of equilibrium properties. 
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I. INTRODUCTION 



One modern approach to therapeutics is to identify and to block the molecular interac- 
tions responsible for disease. Within this context, combinatorial peptide screening experi- 
ments play an important role in the discovery of inhibitors. Libraries of peptides are used 
to discover the crucial molecular structure, or pharmacophore, that can block effectively 
the aberrant binding event. Cyclic peptides are preferred for this purpose since they display 
higher binding affinities due to their reduced conformational entropy in solution.''' The cyclic 
scaffolds and templates derived from combinatorial cyclic peptide studies have been used to 
assemble various spatially-defined functional groups, and highly active analogs have been 

synthesized in this way.'"'"''' A classic example is the blocking of the molecular event respon- 
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sible for blood platelet aggregation by peptides of the form CRGDxxxC, CxxxRGDC, and 

CxxxKGDC.'' Organic analogs of these RGD peptides have been synthesized, and several 
companies are pursuing therapeutic applications in clinical trials. 

The ability of a peptide molecule to bind selectively to a receptor depends on its structural 
and conformational properties, which are in turn determined by its constituent amino acids. 
As the only natural secondary amino acid, proline plays a particular role in peptide and 
protein structural biology as a turn-promoting unit.^ Although most amide bonds in native 
peptides and proteins are in the trans state, a significant fraction (6%) of the X-Proline 
amide bonds are in the cis state.' In some proteins' and cyclic peptides,' prolyl cis/trans 
isomerization has been detected in solution. This isomerization event is frequently a rate- 
limiting step in protein folding. The free energy barrier for the cis/trans isomerization of 
typical X-Pro amide bonds is about 19 Kcal/mol.' While isomerization of non-prolyl amide 
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bonds is rare, the energy barrier to isomerization is similar.- 

Conventional molecular simulations cannot follow the natural cis /trans isomerization 
dynamics, since the typical equilibration time is in the range of 10 to 100 s.- Monte Carlo 
simulations have been successful only in the case where isomerization is experimentally 
known to occur. In such cases, an isomerization reaction coordinate can be defined, and 
special techniques such as umbrella sampling can be applied. Conventional molecular dy- 
namics for cyclic peptides, even in the absence of cis /trans isomerization, turns out to be 
non-ergodic as well, failing to sample the multiple solution conformations during the ac- 
cessible simulation time. The difficulties arise from the cyclic constraints that, along with 
the intrinsic high energy barriers, isolate the accessible conformations to several separated 
regions in phase space. Monte Carlo methods, however, do not require the system to fol- 
low the natural trajectory. Therefore, larger and unphysical moves can be performed to 
overcome these energy barriers. A numerical peptide rebridging method, inspired by the 
alkane rebridging method'' and the configurational bias Monte Carlo (CBMC) method,'"''' 
has been successfully appUed to the simulation of non- proline- containing, cyclic peptides.'' 
Parallel tempering was shown to be a key factor in the efficiency of equilibration. Proline 
poses some geometrical complexity in the rebridging approach, however, and has not yet 
been treated. 

Here we present a new Monte Carlo method, the analytical rebridging scheme, that is 
suitable for equilibration of proline-containing, cyclic peptides. Our analytical method was 
inspired by the solution for a related inverse kinematics problem in robotic control.^ The 
method can accommodate any rigid unit geometry. The rebridging method is not restricted 
to peptides and can be readily applied to other molecules. The cis /trans isomerization is 
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naturally incorporated in the method by allowing for two states of the amide bond. As an 
added benefit, we find the analytical approach to be at least ten times more computationally 
efficient than the previous numerical method, even in the simplest, non-proline-containing 
cases. The analytical rebridging method and other components of our simulation method- 
ology are described in Section |T[ We apply our method to five cyclic peptides that have 



previously been characterized by NMR methods in Section |T|. We show that our method 
can effectively equilibrate these molecules, yielding conformations consistent with the NMR 
analyses. We discuss the results in Section |^ and conclude in Section 0. 



II. SIMULATION METHODS 

A. Analytical Peptide Rebridging 

Our equilibration scheme involves three types of moves: rebridging moves, semi-look- 
ahead (SLA) moves for side chains, and swapping moves. The rebridging move is described 
below and in Sec. |llij|. We describe the SLA move in Sec. fTO and the swapping move in 



Sec. [rm 



In our system, bond lengths and bond angles are kept at their equilibrium value. We 
focus on sampling the biologically-relevant, torsional degrees of freedom. With this simpli- 
fication, a molecule is comprised of a set of rigid units. ^ The rebridging scheme can easily 
be generalized to accommodate flexible bond angles and bond lengths. 

A peptide rebridging move causes a local conformational change within the molecule, 
leaving the rest of the molecule fixed. Consider the segment of a peptide backbone shown in 
Fig. 0. The angles 0o and 07 are rotated, causing the rigid units between and 6 to change. 
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The range of rotation is within iA^max- These two rotations break the connectivity of the 
molecule. We then find all the solutions that re-insert the backbone units in a valid way 
between rigid units 1 to 6. The rebridging move is based upon enforcing six geometrical 
constraints, and this is why we choose to modify six rigid units. Modification of more than 
six units in a single move is possible by the rotation of more than two angles. Such a move 
is likely, however, to lead to an infeasible geometry most of the time and so to result in 
ineffective equilibration. Therefore, we choose to rotate 0o a^nd 07 only. 

Our peptide rebridging scheme features an analytical solution of the geometrical problem 
arising from the reconnection of the backbone units, a problem previously solved in a numer- 
ical way.^'' The side chains are rigidly rotated for each of the solutions. For rebridging moves, 
the solutions for both the new and the old configurations are needed so as to satisfy detailed 
balance. The analytical solution involves the reduction of twenty linear equations to an 8x8 
determinant equation of one torsional angle. The determinant equation is equivalent to a 
polynomial of degree sixteen. Therefore, the maximum number of new geometrical solutions 
is strictly limited to 16, a bound that is obeyed in previous simulations using the numerical 
rebridging method.'' The determinant equation is reformulated as an eigenvalue problem^ 
and solved using the QR algorithm.^ The details of the analytical rebridging method are 
described in Appendix 0. 

Following the "with Jacobian" (WJ) biasing,'' one of the solutions is picked with a prob- 
ability proportional to J^"^-* exp(— /3U*^°^), the product of the Jacobian and the Boltzmann 
factor of the solution. The Jacobian J accounts for the correction to the non-uniform dis- 
tribution generated by rebridging moves, and can be expressed in several ways. Here we 
present one form, involving the determinant of a 5 x 5 matrix B: 
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J 



U6 • 63 

det IBI 



Bij = [uj X (re - r^)],, if i < 3, 



[uj X U6],_3, if i = 4, 5. 



(1) 



Here Uj is the unit vector about which the torsional angle 0i is measured, is the position 
of the atom in unit i — 1 that is bonded by a sigma bond to unit i, and 63 is a unit vector 
along the laboratory 2-axis. The subscript outside the brackets refers to the component of 
the bracketed vector in the laboratory frame. The Jacobian can also be written in terms 
of a determinant of a 4 x 4 matrix, although the definition of the components of the 4x4 
matrix is more involved^ The attempted move is accepted with the probability 



where W'-'^^ and W^"-* are the normalization (Rosenbluth) factors for the new and old solu- 
tions, respectively.' 



For each amide bond that we wish to allow to isomerize, we assign two discrete states 



values of cij = 0° or u; = 180° in the cis and trans conformations, respectively. The partition 
function includes a summation over both states. Because of the sum over states, solutions 
corresponding to all the possible cis/trans states between unit and unit 6 are included 
in the calculation of the Rosenbluth factor. The same approach can be applied to both 
non-prolyl and prolyl amide bonds. 




(2) 



B. The cis/trans Isomerization 



to the corresponding rigid unit. As shown in Fig. |^, the amide unit takes the torsional 
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C. Side Groups 



The chemical functionahty of peptides hes mostly in the side chains. Efficient equili- 
bration for the side chains, therefore, is very important. The semi-look-ahead (SLA) move, 
based on CBMC methods, has been shown to equilibrate effectively long and bulky chains.'' 
We use the SLA method to equilibrate side chains and end groups in this work. A SLA 
move proceeds by regrowing a randomly selected side chain, unit by unit, beginning from 
the bond that connects the backbone to the side chain. The reverse move is performed so 
as to satisfy detailed balance. The Jacobian for each solution is unity. The attempted move 
is accepted with the probability 



where W*^"^ and W^"-* are the normalization (Rosenbluth) factors for the new and old ge- 
ometries, respectively.^ 



Parallel tempering was first proposed for the study of glassy systems with large free 
energy barriers.' It has since been successfully applied to a variety of systems. This 
method achieves rigorously correct canonical sampling, and it significantly reduces the equi- 
libration time in a simulation. Instead of a single system, we consider in parallel tempering a 
larger ensemble with n systems, each equilibrated at a distinct temperature T^, i — 1, . . . , n. 
The system with the lowest temperature is the one of our interest; the higher temperature 
systems are added to aid in the equilibration of the system of interest. In addition to the 
normal Monte Carlo moves performed in each system, swapping moves are proposed that 




(3) 



D. Parallel Tempering 
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exchange the configurations between two systems i and j — i + l<i<n. A swapping 
move is accepted with the probabihty 

acc[(i, j) ^ {j, i)] = min[l,exp(-A/3AU)] , (4) 

where A/3 and AU are the difference of the reciprocal temperatures and energies, respectively. 
The higher temperature systems are included solely to help the lowest temperature system to 
escape from local energy minima via the swapping moves. To achieve efficient sampling, the 
highest temperature should be such that no significant free energy barriers are observed. So 
that the swapping moves are accepted with a reasonable probability, the energy histograms 
of systems adjacent in the temperature ladder should overlap. 

III. COMPARISON WITH EXPERIMENTAL RESULTS 

We perform simulations on five distinct cyclic peptides that have been previously charac- 
terized by NMR methods. We focus on the backbone structure of proline-containing, cyclic 
peptides that were observed experimentally to undergo cis /trans isomerization. Molecular 
interactions are described by the AMBER force field with exphcit atoms.- Aqueous solvent 
effects are estimated by simple dielectric theory.'' Five or six systems were used in the par- 
allel tempering for each simulation, with the highest temperatures ranging from 10^ K to 
10^ K. The lowest temperature system in each case is 298 K. For the first three peptides, 
cis /trans isomerization is allowed in prolyl amide bonds only. For the last two peptides, 
isomerization of all amide bonds is allowed. The simulations take 5-8 CPU hours for the 
first three peptides and 15-20 CPU hours for the last two peptides. All the simulations 
were performed on an Intel Pentium II 450 MHz Linux workstation. Rapid equilibration is 
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achieved for all peptides. The results for each of the peptides are presented below: 

1. c(Pro-Phe-D-Trp-Lys-Thr-Phe). This analog of somatostatin displays high activity 
in inhibiting the release of a growth hormone.''''' Analysis of the NMR spectrum in 
D2O solution indicated a unique backbone conformation.' The prolyl amide bond at 
Phe-Pro adopted a cis conformation. In our simulation, we find an essentially unique 
conformation, possessing the same amide bond cis /trans sequence. A representative 
conformation for this peptide is shown in Fig. ^. 

2. c(Phe-Phe-Aib-Leu-Pro). This pentapeptide contains the Pro-Phe-Phe sequence that 
has been proposed to be responsible for the cytoprotective ability of antamanide and 
cyclolinopeptides.''''' NMR analysis indicated that the peptide is conformationally non- 
homogeneous at room temperature. Two predominant cis/trans isomers for the Leu- 
Pro amide bond were identified in acetonitrile at 240 K.' Our simulation led to two 
inter-converting conformers in the simulation, as shown in Fig. ^. The cis and trans 
conformers occur with probability 58% and 42%, respectively. 

3. c(Gly-Glu(OBzl)-Pro-Phe-Leu-Pro). This cyclic hexapeptide was synthesized for use 
as a possible chiral site for enantiomeric separation. '' NMR studies in dimethyl sulfoxide 
(DMSO) reported two isomers, one having two cis prolyl bonds, and the other having 
all-trans bonds. We find only the 2- cis conformer in simulation, as shown in Fig. |^. 
All the torsional angles fluctuate around mean values, except that the amide group 
between Pro and Gly flips between two opposite orientations. The all-trans conformer 
was found at higher temperatures. 

4. c(Pro-Ala-Pro-Ala-Ala). This cyclic peptide has been designed to serve as a rigid 
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structural template. An unique conformation with two cis prolyl amide bonds in 
DMSO solution was found, according to NMR analysis/ The backbone consists of two 
intertwined type-VIb jS turns, centered about the two prolyl amide bonds respectively. 
We find the same unique conformation. The torsional angles are close to the values 
derived from restrained molecular dynamics.' Figure |^ depicts the geometry of this 
peptide. 

5. Tentoxin, c(MeAla-Leu-MePhe[(Z) A]-Gly). This tetrapeptide selectively induces 
chlorosis in the seedlings of plants. Although tentoxin lacks proline, its two methylated 
amide bonds were found to adopt the cis conformation in a nearly saturated aqueous 
solution.^ The other two non-methylated amide bonds adopt the trans conformation. 
The observation of this cis-trans-cis-trans sequence of the backbone, along with other 
experimental data, led to a proposed boat-like conformation, with the two cis bonds 
located on the same side of the mean plane. In our simulation, we find the same 
amide bond sequence and boat-like conformation that was found in the experimental 
structure. The conformation is shown in Fig. ^ All the carbonyl groups lie in the 
same side of the mean plane, which implies that we have found the the third major 
(conformer C, 8% abundance) of the four conformers found at 268 K in ref. ?. The 
four conformers differ only in the orientation of the two non-methylated amide groups. 
The other three conformers were found at higher temperatures in the simulation. 
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IV. DISCUSSION 



By analyzing the energy trajectories and conformational data, we find that all the pep- 
tides are effectively equilibrated with relatively few Monte Carlo steps. For example, we 
find that the cis/trans equilibrium for peptide 2 was attained within the first 10% of the 
simulation time. Parallel tempering is crucial for this inter-conversion, since the peptide es- 
sentially does not undergo cis/trans isomerization in a single, room-temperature, canonical 
simulation. 

For peptide 1, the NMR-based conformational study indicated a type-ir/3 turn in the 
Phe-D-Trp-Lys-Thr region. This turn is characterized by the hydrogen bond between the 
C=0 of Phe and the N-H of Thr. In our simulation. The C=0 of Phe and the N-H of Thr 
are close to each other, but are not in precise alignment. Such disorder is expected at finite 
temperature. 

Bond angle infiexibility may influence the predicted equilibrium properties of these 
highly-strained molecules. We have investigated the dependence on bond angles by chang- 
ing the angle between 0^-0 and N-Cq, in the amide bond on peptide 2 from 6° to 0°. The 
predicted cis/trans equilibrium shifts from 42% trans to 14% trans, which is a non-trivial, 
although energetically small, cff'ect. We suspect, therefore, that the absence of the all-^rans 
conformer of peptide 3 in our simulation may be due partly to inflexibility of the bond an- 
gles. We used a rigid proline ring, with a 0Pro = —75°. This constraint suppresses the small 
fluctuations that occur in the proline ring. In exceptionally-constrained, cyclic peptides, ^Pro 
can take on other values. For example, ^Pro — —50° has been observed for some unusual 
prolines in the trans state. ^ Fluctuations in the proline ring may, therefore, be important in 
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some highly-strained systems. 

Although experimental conformational analysis of peptide 3 in water was not performed, 
NMR analysis in CHCI3 indicated a minor conformer in addition to the two major ones 
found in DMSO solution.' It is clear, then, that solvent effects play a role in the equilibrium 
structure of peptide 3. Note that the dielectric constants are ecHcis = 4.8,- £dmso = 45,- 
and ShzO = 78.^ The low-dielectric CHCI3 solution favored the all-^rans conformer. Water, 
with a high dielectric constant, may favor the 2-cis conformer, which is what we observe in 
simulation. Indeed, upon reducing the dielectric constant in our implicit solvent model, we 
find a small amount of the all-trans state. However, the detailed structure of the solvent 
molecules around the peptide is likely to be important, and a more accurate description of 
solvent is likely necessary to account fully for the solvent effects. 

Tentoxin (peptide 5) in water at 268 K was experimentally found to aggregate in a 
way that suggested micellar organization. The critical micelle concentration (CMC) was 
estimated to be roughly 35 /iMJ Conformations observed at concentrations above the CMC 
may differ from that of the dilute-limit, monomeric form. In fact, of all four conformers 
found in ref. ?, only conformer C yielded chemical shifts for the 5 and 7 protons of Leucine 
close to that expected for a monomeric form. The other two major conformers, A and B, 
displayed strong shielding effects. The shielding constants for the minor conformer D were 
not reported. These shielding effects were explained by an aggregated structure of these two 
conformers. Since the concentration of conformer C (250/iM) is still well above the critical 
micelle concentration, it is not entirely clear why the chemical shifts for conformer C were 
relatively unaffected by potential aggregation.' Nonetheless, the experimental data suggests 
that conformer C may be either in monomeric form or in an environment similar to aqueous 
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solution. This gives one possible explanation for the absence of conformers A and B in our 
simulation, since these conformers are certainly not in monomeric form. 

V. CONCLUSION 

We presented a new method, the analytical rebridging scheme, for the simulation of 
chemically-diverse, chain-like molecules. The method naturally accommodates cis/trans 
isomerizations. Our rebridging scheme, combined with parallel tempering and biased Monte 
Carlo, is very successful at equilibration of proline-containing, cyclic peptides. Our method 
is not limited to cyclic peptides and can simulate any chain-like molecule. We compared our 
simulations with experimental data on five cyclic peptides and found the predicted confor- 
mations to be reasonably accurate. We were able to sample multiple, relevant conformations 
separated by high energy barriers, a feat not possible with conventional molecular dynamics 
or Monte Carlo. The numerical quality of our predictions, while not limited by sampling 
issues, may be limited by our choice of a simple forcefield. Nonetheless, our method can be 
easily extended to accommodate flexible bond angles and bond lengths. In addition, solvent 
effects may be represented more accurately by better implicit solvent models.''^ 

The methods described here are powerful enough and general enough to influence the 
preferred approach to simulating biological systems. For example, our method should be a 
valuable tool for the fitting of new potential parameters for biological systems. New sim- 
ulation methods for long alkanes have made possible the optimization of force fields that 
significantly reduce the discrepancies between simulation and experiment.'*'^ The same ap- 
proach should lead to improved forcefields for biological systems. In addition, we expect that 
our peptide rebridging scheme, combined with parallel tempering, should replace high tem- 
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perature molecular dynamics and is readily suitable for use in NMR-based conformational 
analyses of biomolecules. 
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APPENDIX A: ANALYTICAL METHOD FOR SOLUTION OF THE 

REBRIDGING PROBLEM 



Our analytical method was inspired by the inverse kinematics problem of six-revolute 
manipulators, which is important for automatic control of robotic arms. This problem 
has been proved to have at most sixteen solutions in the general case.'' Lee and Liang 
reduced the problem to a polynomial of one variable of degree sixteen.'''' The polynomial 
is derived by equating the determinant of an 8x8 matrix to zero, each element being a 
quadratic polynomial. A different closed form was obtained later by Raghavan and Roth.'' 
An excellent review of this subject is given by Manocha and Canny.'' In the following, we 
apply the symbohc formulation of Lee and Liang to reduce our peptide rebridging problem 
to an eigenvalue problem. 

The 6-revolute inverse kinematics problem can be formulated as an equivalent closed, 7- 
revolute mechanism, as shown in Fig. ^ The closed loop consists of seven joints with offsets 
ui, U2, . . . , U7 and of seven links with vectors ai, a2, . . . , ay. For a given backbone fragment 
to be rebridged, the corresponding closed loop is uniquely determined by calculating the link- 
joint intersections. We first draw the seven joint lines that are parallel to the incoming sigma 
bonds of the seven rigid units. The links are defined as the shortest vectors that connect 
consecutive joint lines. Therefore, each link is perpendicular to the two adjacent joint 
vectors. We denote the unit axes of the joint Uj and the link aj by and a^, respectively. 
The joint rotation angles 0i, . . . , 07 are the torsional angles measured around the joints; 
they are equal to the corresponding biological torsional angles plus constant offsets. The 
joint rotation angle 07 and other parameters are determined by the given backbone geometry. 
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and we need to calculate only the six unknown joint rotation angles. 

The idea is to find an over-constrained set of twenty equations that are linear in xq = 
tan 06/2 and the sines of cosines of 0i, 02, 04, and 05. These equations are obtained by 
equating scalar and vector products of the loop axes in both directions of the loop. For 
instance, we equate U3 • U6(04, 05) = Ue ■ U3(0i, 02) in the first equation (see below). These 
equations were first derived using a recursive notation.'' We define the chain vector R"' 
as the vector summation of the consecutive joints and links from unit a to unit (3. The 
summation always goes in the direction of units 1,2, . . . , 7. The indices a and (3 can take 
one of two forms: i denotes either starting from or ending at Uj, and i' denotes either starting 
from or ending at aj. The index (3 can be less than a, and this indicates wrapping around 
the closed loop. For example, 

R^'' ^ = ae + U7 + ay + ui + ai + U2 

R^'' * = . (Al) 

The twenty equations are listed below. The left-hand side of each equation is a linear 
function of cos 01, sin0i, cos 02, sin 02, and Xq. The right-hand side of each equation is a 
linear function of cos 04, sin 04, cos 05, sin 05, and 

U3 • Ue = Ue • U3 

U3 -Ue Xq^Uq- U3 Xq 

R^'' ^' • U6 X U3 = R*^'' ^' • U3 X ue 

R^'' ^' • U6 X U3 X6 = R^'' ^' • U3 X U6 Xq 

''-us^ -R' ^' ■ U3 
R=^' ^' -113X6 = -R6' 2' . 
16 



R^^' • U6 = -R'^' • U6 



6, 2' 



R^' ^' • U6 a;6 = -R^' ^' • Ue 



R3' 5' . r3- 5' ^ j^6, 2' . j^6, 2' 



R^' 5' • R=^' Xe = R^' • R^' 



X6 



Us • as - Us • Ue X as = -ae • U3 Xe - Ue x ae • U3 



U3 • Ue X as Xe + U3 • as = -Ue x ae • U3 xe + ae • U3 



-R^' ^ • as Xe + R ' • Ue x as = -R ' • ae xe - R ' • Ue x ae 



R^' ^' ■ Ue X as xe + R'^' • as = R*'' ^' • Ue x ae Xe - R"' ■ ae 



.3, 5' 



6, 2' 



.6, 2' 



R^'' ^ • as X U3 Xe - R^'' ^' • (ue x as) x U3 

= -(R^' ^' • U3 X ae - Ue X ae • U3) xe - (R^'' ^' • U3 x (ue x ae) + |ue|ae 



R 



3', 5' 



(ue X as) X U3 Xe + R^'' ^ • as x U3 



= -(R^'' ^' • U3 X (ue x ae) + juejae • U3) xq + (R^' ^' • U3 x ae - ue x ae 
1/2(R='' . ''){u, ■ ue) - (R=^' ■ U3)(R'' ■ ue) 
= 1/2(R^' . R^' ^')(ue • U3) - (R^' ■ Ue)(R^' ■ U3) 

1/2(R3' 5' . R3' 5')(u3 . Ue) - (R3' 5' . U3)(R3. 5' . 



1/2(R«' . R«' ^')(ue • U3) - (R^- ■ U6)(R'' ■ U3) 



.6, 2' 



Xe 



Xe 



1/2(R^' . ^')(U3 ■ as) - (R^' ■ k,){R'' ■ U3) 

- [l/2(R^' 5' . R=^' ^')(U3 • Ue X as) - (R^' ■ U3)(R=^' ^' • Ue x as) 
= - [i/2(R6' 2' • R6' 2')(u3 . ae) - (R^' ■ ae)(R^' ■ U3)] ^e 

- [i/2(R6' 2' . j^6, 2')(^^ . Ue X ae) - (R^' ■ U3)(R^' ■ Ue 
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X 



1/2(R3- 5' . r3. 5')(u3 . ug X as) - (R^- ■ U3)(R'' ■ ue x 



3, 5' 



>3, 5' 



+ 



1/2(R6' 2' . r6. 2') . X as) - (R^- 2' ■ U3)(R'' ■ ue x ag 
1/2(R6' 2' ■ R^' ^')(U3 ■ ae) - (R^' ■ k,){-R'' ■ us) 

This set of equations can be put into matrix form as 



(A2) 



AB = Px6 + Q , 



(A3) 



where A is a 20 x 16 constant matrix, and B is a 16 x 1 matrix with the following variables 
as its elements: 



B^ — [C4C5X6, S4C5XQ, c^Xq, C4S5X6, S4S5X6, S5X6, C4X6, S4X6, 



C4C5, S4C5, C5, C4S5, S4S5, S5, C4, S4J . 



(A4) 



Here q = cos0j and Si = sin^j. The 16 x 1 matrices P and Q have elements that are linear 
in the sines and cosines of 0i and 02- Sixteen equations are chosen from the twenty equations 
in eq. (|A2|) and used to express B as functions of 0i, 02, and Xq. This is accomplished by 



defining As as the corresponding 16 x 16 sub-matrix of A and multiplying both sides of 
these sixteen equations by the inverse of As to obtain 



B = As-^PsX6 + As-^Q 



-1/ 



(A5) 



Here Ps and Qs are the corresponding sub- matrices of P and Q. We find that special 
geometries, such as vanishing lengths of links, may render some choices for Ag singular. The 
linear dependence is identified by singular value decomposition- and avoided by choosing 
16 linearly-independent equations. The expression for B is inserted into the other four 



unused equations to yield four equations that are linear functions of xq and the sines and 
cosines of 0i and 02- Replacing the sines and cosines of 0i and 02 with xi = tan0i/2 and 
X2 = tan 02/2, these four equations can be expressed as 



{diX2^ + eiX2 + fi)x6 + {giX2^ + hiX2 + Pi) = 0, i = 1, . . . , 4, 



(A6) 

where rfj, e,, . . . , Pi are quadratic in Xi. Multiplying eq. ( |A6|) by X2 gives four additional 
equations, which together with eq. (|A6|) can be used to set up a linear set of equations: 
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(A7) 



Here is an 8 x 1 null vector. Solutions are found by equating to zero the corresponding 
determinant of the 8x8 matrix, which can be expanded to a polynomial of degree sixteen in 
Xi- The determinant equation is best solved by reformulating it as an eigenvalue problem. 
This is achieved by rewriting eq. ([A7|) as 



Ao + AiXi + A2a;i^ 







(A8) 



where Aq, Ai, and A2 are 8x8 numerical matrices. Note that in the crankshaft case, Aq, 
Ai, and A2 are identically zero, and so any solution for xi is possible.^ The roots of the 
determinant equation are the eigenvalues of the matrix^ 



I 

-A2"^Ao -A2"^Ai 



(A9) 



where and I are 8x8 null and identity matrices, respectively. The matrix is first reduced 
to an upper Hessenberg matrix, and then the eigenvalues are found by the QR algorithm"" 
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to obtain xi. Occasionally the matrix A2 may be almost singular. This occurs when one of 
the roots is 0i ~ ivr. This singularity can be avoided by the transformation 

tiXi + t2 / A 1 n^ 

^1 = , _ . , , (AlO) 

where ti, ^2, ^3, and are random numbers uniformly distributed in [-1,1]. Solutions for 



Xi are substituted into eq. ( [A^ ) to calculate X2 and xq. These values are substituted back 
to eq. (|A5|) to calculate ^4 and ^5. The angle 03 is not needed for our purpose; it can be 
determined by calculating the positions of a2 and as from the other torsional angles. 



FIGURES 

FIG. 1. A backbone segment selected to be rebridged. Only the backbone atoms are shown. A 
change of the driver angles 0o and ^7 breaks the connectivity. The dotted area represents the region 
in which the positions of the backbone atoms must be restored. The thick solid lines represent pi 
bonds or rigid molecular fragments within which no rotation is possible. 

FIG. 2. The (a)c«s and {h)trans conformation of the amide bond. Only backbone atoms are 
depicted. The Rx atom is hydrogen for a normal, non-prolyl amide bond or carbon for a prolyl or 
methylated amide bond. 

FIG. 3. The conformation of c(Pro-Phe-D-Trp-Lys-Thr-Phe). The colors green, red, and blue 
denote carbon, oxygen, and nitrogen atoms, respectively. The hydrogen atoms are omitted. 

FIG. 4. The (a)cis and (h)trans conformations of c(Phe-Phe-Aib-Leu-Pro). 

FIG. 5. The conformation of c(Gly-Glu(OBzl)-Pro-Phe-Leu-Pro). 

FIG. 6. The conformation of c(Pro- Ala-Pro- Ala- Ala). 

FIG. 7. The conformation of c(MeAla-Leu-MePhe[(Z) A]-Gly). 

FIG. 8. The geometry of the closed, 7-revolute mechanism, consisting of 7 joints and 7 links. 
The joints are represents by Ui, U2, . . . , uy. The links are represented by ai, a.2, ■ ■ ■ , ay. Each 
link is perpendicular in three dimensions to the two adjacent joints. The unit axes of the joints 
and links are defined as Uj and a^, respectively. 




Figure Wu and Deem, "Analytical Rebridging Monte Carlo. 
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MePhe[(Z)A] 
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